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Abstract 



In this paper we consider a discrete scale invariant (DSI) process {X(t), t e R + } with 
scale / > 1. We consider to have some fix number of observations in every scale, say T, 
and to get our samples at discrete points a k , k £ W where a is obtained by the equality 
I = a T and W = {0, 1, . . .}. So we provide a discrete time scale invariant (DT-SI) pro- 
cess X(-) with parameter space {a k ,k g W}. We find the spectral representation of the 
covariance function of such DT-SI process. By providing harmonic like representation of 
multi-dimensional self-similar processes, spectral density function of them are presented. 
We assume that the process {X(t),t e R + } is also Markov in the wide sense and pro- 
vide a discrete time scale invariant Markov (DT-SIM) process with the above scheme 
of sampling. We present an example of DT-SIM process, simple Brownian motion, by 
the above sampling scheme and verify our results. Finally we find the spectral density 
matrix of such DT-SIM process and show that its associated T-dimensional self-similar 
Markov process is fully specified by {Rf (1), Rf (0), j = 0, 1, . . . , T - 1} where Rf(r) is 
the covariance function of jth and (j + r)th observations of the process. 
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1 Introduction 

The concept of stationarity and self-similarity are used as a fundamental property to han- 
dle many natural phenomena. Lamperti transformation defines a one to one correspondence 
between stationary and self-similar processes. A function is scale invariant if it is identical 
to any of its rescaled version, up to some suitable renormalization in amplitude. Discrete 
scale invariance (DSI) process can be denned as the Lamperti transform of periodically corre- 
lated (PC) process [2]. Many critical systems, like statistical physics, textures in geophysics, 
network traffic and image processing can be interpreted by these processes [Tj. Flandrin et. 
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al. introduced a multiplicative spectral representation of DSI processes based on the Merlin 
transform and presented preliminary remarks about estimation issues [2], [6]. 

As the Fourier transform is known as a suited representation for stationarity, but not for 
self-similarity. A harmonic like representation of self-similar process is introduced by using 
Merlin transform [6]. 

Covariance function and spectral density of discrete time periodically correlated Markov 
process have been studied and characterized in [llj . Markov processes have been the center of 
extensive research activities and wide sense Markov processes are studied before the general 
theory. In some texts, these processes are defined in the case of transition probabilities of 
a Markov process. Various classes of wide sense Markov processes are like jump processes, 
diffusion processes and processes with a discrete interference of chance [7J. A process which 
is Markov and self-similar, is called self-similar Markov process. These processes are involved 
in various parts of probability theory, such as branching processes and fragmentation theory 

®- 

In this paper, we consider a DSI process with some scale / > 1, and we get our samples 
at points a k , where k € W, I = a T , W = {0, 1, . . .} and T is the number of samples in each 
scale. By such sampling we provide a discrete time scale invariant process in the wide sense 
and find the spectral representation of the covariance function of such process. 

This paper is organized as follows. In section 2, we present stationary and self-similar 
processes by shift and renormalized dilation operators. Then we provide a suitable platform 
for our study of discrete time self-similar (DT-SS) and discrete time scale invariant (DT-SI) 
processes by introducing quasi Lamperti transformation. Harmonizable representation of 
these processes are expressed in this section too. Also by using the spectral density matrix of 
PC processes, the spectral representation of the covariance function of DT-SI processes are 
given. In section 3, a harmonic like representation of multi-dimensional self-similar processes 
and spectral density function of them are obtained. As an example we introduce a process 
called simple Brownian motion which is DSI and Markov too. Finally a discrete time scale 
invariant Markov (DT-SIM) process with the above scheme of sampling is considered in 
section 3 and the spectral density matrix of such process and its associated T-dimensional 
self-similar Markov process are characterized. 



2 Theoretical framework 

In this section, by using renormalized dilation operator, we define discrete time self-similar 
and discrete time scale invariant processes. The quasi Lamperti transformation and its prop- 
erties are introduced. We also present the harmonizable representation of stationary and 
harmonic like representation of self-similar processes. The spectral density of PC processes 
and the spectral representation of the covariance function of DT-SI processes are given. 
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2.1 Stationary and self-similar processes 

Definition 2.1 Given r G R, the shift operator S T operates on process {Y(t),t G R} ac- 
cording to 

(S T Y)(t):=Y(t + r). (2.1) 
A process {Y(t),t G R} is said to be stationary, if for any t, r G R 

{(S T Y)(t)} I {Y(t)} (2.2) 

where = is the equality of all finite- dimensional distributions. 

If (2.2) holds for some r G R ; i/ie process is said to be periodically correlated. The smallest 
of such t is called period of the process. 

Definition 2.2 Given some numbers H > and A > 0, the renormalized dilation operator 
T^H,x operates on process {X(t),t G R + } according to 

(V HiX X)(t):=\- H X(\t). (2.3) 

A process {X(t),t G R + } is said to be self-similar of index H , if for any A > 

{{V H:X X){t)} ± {X(t)}. (2.4) 

The process is said to be DSI of index H and scaling factor Xq > or (H,Ao)-DSI, if (2.4) 
holds for A = Ao • 

As an intuition, self-similarity refers to an invariance with respect to any dilation factor. 
However, this may be a too strong requirement for capturing in situations that scaling prop- 
erties are only observed for some preferred dilation factors. 

Definition 2.3 A process {X(k),k G T} is called discrete time self-similar (DT-SS) process 
with parameter space T, where T is any subset of distinct points of positive real numbers, if 
for any k\ , /c2 £ T 

{X(k 2 )} I (h)H {x{ki)} (2 5) 

The process X(-) is called discrete time scale invariance (DT-SI) with scale I > and param- 
eter space f, if for any k\,k 2 = lk\ G T, (2.5) holds. 

Remark 2.1 // the process {X(t),t G R + } is DSI with scale I = a T for fixed T G N and 
a > I, then by sampling of the process at points a k ,k G W where W = {0, 1, . . .}, we have 
X(-) as a DT-SI process with parameter space T = {a k ;k G W} and scale I = a T . If we 
consider sampling of X(-) at points a nT+k , n G W,/or fixed k = 0, 1, . . . , T — 1, then X(-) is 
a DT-SS process with parameter space T = {a nT+h ;n G W}. 
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Yazici et.al. [13] and [13] introduced wide sense self-similar processes as the following 
definition, which can be obtained by applying the Lamperti transformation Ch to the class 
of wide-sense stationary processes. This class encompasses all strictly self-similar processes 
with finite variance, including Gaussian processes such as fractional Brownian motion but no 
other alpha-stable processes. 

Definition 2.4 A random process {X(t),t G R + } is said to be wide sense self-similar with 
index H, for some H > if the following properties are satisfied for each a > 

(t) E[X 2 (t)] < oo, 

(ii) E[X(at)] = a H E[X{t)], 

(Hi) E[X(at 1 )X(at 2 )] = a 2H E[X(ti)X(t 2 )) . 

This process is called wide sense DSI of index H and scaling factor ao > 0, if the above 
conditions hold for some a = a$. 

Definition 2.5 A random process {X(k),k E T} is called DT-SS in the wide sense with 
index H > and with parameter space T, where T is any subset of distinct points of positive 
real numbers, if for all k,k\ € T and all a > 0, where ak, ak\ £ T : 

(i) E[X 2 (k)] < oo, 

(ii) E[X(ak)\ = a H E[X(k)}, 

(Hi) E[X(ak)X(ak 1 )] = a 2H E[X(k)X(k 1 )\. 

If the above conditions hold for some fixed a = a®, then the process is called DT-SI in the 
wide sense with scale ao- 

Remark 2.2 Let {X(t),t 6 R + } in Remark 2.1 be DSI in the wide sense. Then X(-) with 
parameter space T = {a k ; k E W} for a > 1 is DT-SI in the wide sense, where W = {0, 1, . . .} 
and X(-) with parameter space T = {a nT+k ;n € W} for fixed T G N and a > 1 is DT-SS in 
the wide sense. 

Through this paper we are dealt with wide sense self-similar and wide sense scale invariant 
process, and for simplicity we omit the term "in the wide sense" hereafter. 

2.2 Quasi Lamperti transformation 

We introduce the quasi Lamperti transformation and its properties by followings. 

Definition 2.6 The quasi Lamperti transform with positive index H and a > 1, denoted by 
£-H,a operates on a random process {Y(t),t G R} as 

£ H ,aY(t)=t H Y(log a t) (2.6) 
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and the corresponding inverse quasi Lamperti transform C H a on process {X(t), t G R + } acts 
as 

C H ) a X{t) = a-^a'). (2.7) 

One can easily verify that CH,a£']i~ a X(t) = X(t) and C^ a CH,o¥(t) = Y(t). Note that in the 
above definition, if a = e, we have the usual Lamperti transformation Ch- 

Theorem 2.1 The quasi Lamperti transform guarantees an equivalence between the shift 
operator S\ oga k and the renormalized dilation operator T>jj,k in the sense that, for any k > 

£~H,cPH,k£>H,a = S\og a k- (2-8) 

Proof: 

C- H ] a V Htk C H , a Y{t) = C- H ] a V Htk {t H Y{\og a t)) = C- H ] a {k~ H (kt) H Y(\og a kt)) 
= C H ) a (t H Y(\og a kt)) = cr tH {o?) H Y (log, ka l ) = Y(log a k + t)= S loga k Y(t).D 

Corollary 2.1 If{Y(t),t € R} is stationary process, its quasi Lamperti transform {Cu,aY{t), 
t 6 R + } is self-similar. Conversely if{X(t),t £ R + } is self-similar process, its inverse quasi 
Lamperti transform {C^ a X(t),t € R} is stationary. 

Corollary 2.2 If {X(t),t G R+} is (H,a T )-DSl then Cjf X(t) = Y{t) is PC with period 
T > 0. Conversely if '{Y(t),t G R} is PC with period T then C H , a Y "(t) = X(t) is (H, a T )-DSI. 

Remark 2.3 If X{-) is a DT-SS process with parameter space T = {l n ,n G W} ; then its 
stationary counterpart Y(-) has parameter space T = {nT,n G N}: 

X(l n ) = C H , a Y(l n ) = l nH Y{\og a a nT ) = a nTH Y(nT). 

Also it is clear by the following relation that if X(-) is a DT-SI process with scale I = a T , 
T G N and parameter space T = {a k ,k G W}, then Y(-) is a discrete time periodically 
correlated (DT-PC) process with period T and parameter space T = {n,n G N}: 

Y(n) = C-) a X{n) = a' nH X{a n ). 
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2.3 Harmonizable representation 

A stationary process Y(t), EY(i) = 0, can be represented as 

POO 

Y(t)= / e iut d<p(u) (2.9) 



— oo 



which is called harmonizable representation of the process, and ip(u>) is a right continuous 
orthogonal increment process, see [9]. Also the covariance function can be represented as 



OO fOO 

ituj—isio' 



Ry(t,s)= / / e ltw - %stJ d®(u,u') (2.10) 

J — oo J —oo 

where the spectral measure satisfies 

and d^f(ui) = E[\d(p(u)\ 2 ]. All the spectral mass is located on the diagonal line to = oj' . When 
5>(u;,u/) is absolutely continuous, we have spectral density such that d$(u,uj') = 

cf)(u),U)')du}duj' . A necessary and sufficient condition for this equality to hold, as Loeve's condi- 
tion for harmonizability, is that <3?(u;,u/) must satisfy J J |d<3?(u;, < oo. The corresponding 
notion for processes after a Lamperti transformation introduces a new representation for a 
class of processes deviating from self-similarity, which is called multiplicative harmonizability. 
A self-similar process X{t) has harmonic like representation as an inverse Mellin transform, 
namely an integral of uncorrelated spectral increments d(p(oj) on the Mellin basis pQ. 

X(t) = j t H+iul d<p{uj) (2.12) 

and the process has this property if it verifies as 

R x (t,s)= [ ft H - Hu s H - ib/ d^(u,u/). (2.13) 



The inverse Mellin transformation gives the expression of the spectral function if the corre- 
lation is known as [2] 

= J Jr H ^s- H +^R x (t,s)j^-. (2.14) 
2.4 Spectral density function 

The spectral density of a PC process is introduced by Gladyshev in [8]. If Y(n) is a DT- 
PC process, the spectral density matrix is Hermitian nonnegative definite T x T matrix of 
functions f(u) = [fjk(^)]j,k=o,i,...,T-i, and the covariance function has the representation 

T-l 

R n (r) := Cov{Y(n),Y(n + r)) = £ B k (r)e 2kmn ' T (2.15) 

k=0 
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where 2 ^ 

B k {r)= [ ^ j™f k { U )du. 
Jo 

Also f k (ui) and fjk(w),j, k = 0, 1, . . . , T — 1 are related through 

f jk (co) = l/fc.j- ((w - 2nj)/T) , 0^u<2tt. 

For fc<0, w<0orw> 2tt, the functions are defined by the equality f k+ x(co) = /fc(w) 

and f k (u) = f k (u + 2vr). 

Let {X(t),t G R + } be a zero mean DSI process with scale I. If I < 1, we reduce the time 
scale, so that Z in the new time scale is greater than 1. Our sampling scheme is to get samples 
at points a k , k G W, where by choosing the number of samples in each scale, say T G N, we 
find a by / = a T . Therefore the process under study {X(a n ),n G W} is DT-SI with scale 
I = a T . 

Proposition 2.2 If X(a n ) is DT-SI with scale I = a T , T G N ; then we have the spectral 
representation of the covariance function of the process as 

T-1 

rH(j) ■- Cov(X(a n ),X(a n+T )) = a^ 2n+T)H B k ( T )e 2k7Tin/T (2.16) 

k=0 

WhGre B k (r)= Pe™f k {u)<L, (2.17) 



Jo 

and 

/ifcH = ^/fc-i((w-27Tj)/T) (2.18) 
for j, k = 0, 1, . . . , T - 1 and < u < 2tt. 



Proof: According to (2.6) and Corollary 2.1, for any n,r G W 

R%{t) = E[X{a n )X{a n ^)] = E[C H , a Y( a n )C H , a Y(a n ^)] 

= a {2n+T)H E[Y(n)Y(n + r)] 
where Y(n) is DT-PC process with period T = log a I. Thus by (3.1) 

R H(t) = a( 2n ^ H R n (r) = a^ 2n+ ^ H £ B k {r)e 2kmn / T U 

k=0 
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3 Characterization of the spectrum 



In this section we provide spectral density matrix of multi-dimensional self-similar process 
W(n). By using harmonic like representation of a self-similar process, we characterize the 
spectral density matrix of DT-SI process in subsection 3.1. A discrete time scale invariant 
Markov (DT-SIM) process with a new scheme of sampling is considered and the properties 
of an introduced example is verfied too. The spectral density matrix of such process and its 
associated T-dimensional self-similar Markov process are characterized in subsection 3.2. 



3.1 Spectral representation of multi-dimensional self-similar process 

By Rozanov [12], if = {£, k {t)}k=i....,n be an n-dimensional stationary process, then 



£(t) = J e lM <P(d\) (3.1) 

is its spectral representation, where <f> = {(pk}k=i,...,n an d is the random spectral measure 
associated with the fcth component £ fc of the n-dimensional process £. Let 



Bkr(r) = E[£ k (t + t)?(t)], k,r = l,...,n 

and B{t) = [Bk r (r)]k r =i n be the correlation matrix of £. The components of the correla- 
tion matrix of the process £ can be represented as 

B kr (T) = J e iXr F kr (d\), k,r = l,...,n (3.2) 



where for any Borel set A, F kr (A) = E[ip k (A)ip r (A)] are the complex valued set functions 
which are cr-additive and have bounded variation. For any k,r = 1, . . . , n, if the sets A and 
A' do not intersect, E[<pk(A)(p r (A')] = 0. For any interval A = (Ai,A2) when -Ffc r ({Ai}) = 
Fkr({^2}) = the following relation holds 

1 r °° 

F kr (A) = — V B kr (r)e- lXr d\ (3.3) 

\ I g— «A2T _ g— iXir 

= —B kr (0)[X 2 - Ai] + lim — B kr( T ) : 

0<|t|^T 

in the discrete parameter case, and 

F kr (A) = lim — / : B kr (T)dT 

a^too 27T J_ a —IT 

in the continuous parameter case. 
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Using the above results of Rozanov for multi-dimensional stationary processes and using 
the Lamperti transformation, we present the definition of multi-dimensional self-similar pro- 
cess and obtain the properties of the corresponding multi-dimensional self-similar process by 
the following theorem. 

Definition 3.1 The process U(t) = (U°(t), f7 1 (t), . . . , L r<5 ' _1 (t)) is a q- dimensional discrete 
time self similar process in the wide sense with parameter space T, which consists of finite or 
countable many points o/R + , if the followings are satisfied 

(a) {U 3 (-)} for every j = 0, 1, ■ ■ ■ , q — 1 is DT-SS process with parameter space T. 

(b) U % (-) and U 3 (•) for i,j = 0, 1, • • ■ , q — 1 have self-similar correlation, that is 

Cov (U\ts), U j (tr)) = t 2H Cov (U\s), U j (r)). 
where s,r,ts,tr are in T. 

Theorem 3.1 Let W(a k ) = (W°(a k ), W 1 (a k ), ... , W q ~ l (a k )), fceZ.ol be a discrete 
time q-dimensional self-similar process. Then 

(i) The harmonic like representation ofW J (a k ) is 

/•27T 

W j (a k ) = a kH / e iujk dipjuj). (3.4) 
Jo 

where (pj(ui) is the corresponding spectral measure, that E[d(pj(uj)dip r (uj')} = dDf r (u), j,r = 
0, 1, . . . , q — 1 when uj = uj' and is when oj ^ oo' . We call D^ r {uj) the spectral distribution 
function of the process. 

(ii) The corresponding spectral density matrix of {W(a k ), k G Z} is d H (uj) = [<A-f r (oj)]j, r =Q,...,q-i, 
where 

1 oo 

d f» = ^ E *- nH e-^Qf r (a n ) (3-5) 

n=~ oo 

a > 1 and Qf r (a n ) is the covariance function ofW J (a n ) and W r (l). 

Before proceeding to the proof of the theorem we remind that based on our sampling 
scheme at points a k , k G Z, of continuous DSI process with scale I = a T , a G R, T € N. So 
we consider W(-) at points l n = a nT as the corresponding T-dimensional DT-SS process and 
apply this theorem in Lemma 3.4. 
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Proof of (i): W^(a k ) for j = 0, 1, . . . , q — 1 is DT-SS and its stationary counterpart &(k) 
has spectral representation = J Q 27r e lulk d(pj(ui). Thus by (2.6) 

/•27T 

W"'(a fe ) = C H , a e(a k ) = a kH e{k) = a kH / e^dif^u). 

Jo 

Proof of (ii): The covariance matrix is denoted by Q H (n,r) = [Qf r (a n , a T )]j,r=o,...,g-i 
where its elements have the spectral representation 

Qf r {a m , a T ) = E[W j {a m a T )W r (a m )} = a 2mH E[W j (a T )W^{T)} = a 2mH Qf r (a T ). (3.6) 

Also by (3.4) 

r-2-K r-2-K rlix 

Qf r (a T ) = a TH E[ e iu)T 'd^M / d^ r {u')] = a TH e iu)1 dD? r (oj) (3.7) 
Jo Jo Jo 



where E[d(pj{uj)dip r {oj')] = dD^ r {uj) when oj = u' and is when u ^ uj' . 

The spectral distribution function of the correlation matrix Q H (a k ) = [Q^(a fc )]j, r =o,...,g-i 

is 

D H (du) = [Df r (doj)] jtr= o,..., q -i. 
(3.2)-(3.3), (3.7) and appropriate transformation we have 

oo 

D ?M) = - a- nH e-^Qf r (a n )d\. (3.8) 

•> A n=-oo 

Let A = (uj,uj + duo], then we have the spectral density matrix as d H (lo) = [dj^(a;)] :;if . = o,...,g-i 
where 

D H (dio) 1 00 1 ru+du 



n=— oo 



_L V (- — lim - — W -Hn O?(a n ) 



n=—oo 



1 OO oo 

^ £ ((--)(-He— )a-^Qf r K) = - £ «-^e— Q? r (a»). 



27r ^— ' in 2tt 

n=—oo n=— oo 



Existence of dj^(w) follows from part (i) of the theorem as VF fc (a") is the Lamperti counter- 
part of the stationary process £ fc (n), k = 0, . . . , q — !.□ 
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3.2 Spectral density of DT-SIM process 

Let {X(t),t G R} be a DSI process with scale I and Markov in the wide sense. Using our 
sampling scheme described in this section, we assume I and a to be greater than one. Thus 
{X(a n ),n G W} is a discrete time scale invariant Markov (DT-SIM) process with scale 
I = a T . 

Let R(t\,t2) be some function defined on Tx T and suppose that R{t\,t2) ^ everywhere 
on TxT, where 7~is an interval. Borisov [3] showed that the necessary and sufficient condition 
for R(ti,t2) to be the covariance function of a Gaussian Markov process with time space T 
is that 

R{h,t 2 ) = G(min(ii,i2))ir(max(ti,i 2 )) (3.9) 

where G and K are defined uniquely up to a constant multiple and the ratio G/K is a positive 
nondecreasing function on T ■ 

It should be noted that the Borisov result on Gaussian Markov processes can be easily 
derived in the discrete case for second order Markov processes in the wide sense, by using 
Theorem 8.1 of Doob [5]. 

Here we present a closed formula for the covariance function of the DT-SIM process 
and characterized the covariance matrix of corresponding T-dimensional self-similar Markov 
process by theorems 3.2 and 3.3 [101. 

Theorem 3.2 Let {X(a n ),n G Z} be a DT-SIM process with scale I = a T , a > 1, T £ N, 

then the covariance function 

RH( T ) = E[X(a n+T )X(a n )} (3.10) 
where r G Z, n = 0, 1, . . . , T - 1, R^ +T ( T ) = a2THR n( T ) and Rni T ) + is °f the form 

R%(kT + v) = [/i(a T - 1 )]^(a t ' +n - 1 )[/i(a n - 1 )]- 1 J R^(0) (3.11) 

Rn(-kT + v) = a~ 2kTH Rn +V {(k -l)T + T-v) 
where k £ Z, v = 0, 1, . . . , T — 1, 

r r 

h( a r ) = J] h{cP) = J] Rf(l)/Rf(0), r G Z (3.12) 

3=0 j=0 

and h(a~ 1 ) = 1. 

Proof: Here we present the sketch of the proof. From the Markov property (3.9), for a > 1, 
we have that 

RH(j) = G{a n )K{a n+T ) r G W. (3.13) 
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and 
Thus 



RH(0) = G(a n )K(a n ). 



K{a n+T ) = ^f\K{a n ). 



itf(O)" 



For r = 1, by a recursive substitution in (3.14) one can easily verify that 

n-1 

K(a n ) = K(l) Yl h(a j ) 

3=0 

where h(a?) = Rf(l)/Rf(0). Hence for n = 0, 1, . . . , T - 1, k € W 

fcT+n-l 

K(a kT+n ) = K(l) '[I h(a j ). 

3=0 

As X(-) is DT-SI with scale a T by (3.10) 



( } ^(0) " Rf (o) " * (a } ' 



i e W. 



Therefore using (3.12) we have 



fcT+n-l 

Y\ h{a j ) = [h(a T - 1 )] k h(a n - 1 ). 

3=0 



Consequently for n = 0, 1, . . . , T — 1 

K{a kT+n ) = K(l) 



h(a r - L ) 



n— In 



Let t = kT + v, then it follows from (3.14) and (3.17) that 



R*{kT + v) = 



K{a n+kT+v ) H K(l)[h{a T - 1 )] k h 



-Rn(0) = 



a 



v+n— 1 N 



-^(0) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



= [/i(a T - 1 )] fc /i(a t,+ri - 1 )[/i(a n - 1 )]- 1 ^(0) 

for k = 0, 1, . . ., a > 1 and n, v = 0, 1, . . . , T — 1. Similar to the above proof for r = — kT + u 
we have 

R%(-kT + v) = [^(a T - 1 )]-^K + "- 1 )[^(a"- 1 )]- 1 J R^(0) 
and also note that 

R%(-kT + v) = E[X(a- kT+n+v )X(a n )} = a~ 2kTH E[X{a n+v )X(a kT+n )] 
= a~ 2kTH R^ +v (kT + v) = a- 2kTH R»((k — 1)T + T — v).D 
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Example 3.1 We consider moving of a particle in different environment A±,A2,... based 
on Brownian motion with different rates. Specially we consider this movement by X(t) with 
index H > and scale A > 1 as 

oo 

X(t)=^A B ( H -^/ [A n-l iA »)(t)S(t) 

n=l 

where B(-), /(•) are Brownian motion and indicator function respectively and we call this 
process, simple Brownian motion. 

Let Ai = [1,A), A 2 = [A, A 2 ) and A n = [A n_1 ,A n ) as disjoint sets. The process X(t) is 
DSI and Markov too. For checking these properties, first we find the covariance function of 
it. The covariance function of the process for t £ A n , s £ A m and s ^ t is 

Cov(X(t), X(s)) = A (n+m)(H -^Cov(S(t), B(s)) = \( n + m )( H -\) s ( 3 .1 8 ) 

since as we know Cov (B (t), B(s)) = min{i, s}. Therefore by the condition (3.9), the above 
covariance is the covariance function of a Markov process. Now we verify the DSI property. 
If t is in (A™" 1 , A n ] then At is in (A n , A n+1 ]. Thus for t £ A n+1 and s £ A m+1 we have 



Cov(X(Xt),X(Xs)) = X^ n+m+2 ^ H -^Cov(B(Xt),B(Xs)) = \(n+ m +2)(H-i) Xs 

= X^ x {n +m) {H-^) s = X* H Cov(X(t),X(8)). 

Then X(t) is (H,A)-DSI. 

By sampling of the process X(-) at points a n , n G W, where A = a T , T € N and A > 1, 
we provide a DT-SIM process and investigate the conditions of Theorem 3.2. For j = kT + i 
where i = 0, 1, . . . , T — 2 and k = 0, 1, . . . by (3.18) we have that 

j) - R f {l) - Cov(^(« J+1 ),^(« i )) _ q^+Prg'+J _ 
{a } ~ Rf(0) ~ Cov(X(a3),X(a3)) ~ a ^+i)TW +j ~ 1 

as a j ,a j+1 £ A k+1 and H ' = H - \ . Also for j = kT + T - 1 we have that 

, r = gfO) = Cov(X(^+ 1 ),A-(^)) = q(Wfi'+j = TRI 
[a } Rf(0) Cav{X(Qj),X(a3)) a ^k+2)TW +j a 

as a j £ A k+l and £ A k+2 . Thus for j = kT + i, i = 0, 1, . . . , T - 2 and k = 0, 1, . . . 

kT-\-i k 

h(a kT+l ) = 11 h(a r ) = f[a TH ' = a kTH ' 
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and for j = kT + T - 1 

fcT+T— 1 

/ l ( a ^+ r - i )= n ^)=n« Tif '=« (fc+i)ra '- 

r=0 r=0 

V T-U _ „,TH' 



1 U+n-1 ^T-2 



Finally as h(a 2 1 ) = a 1 

h(a v+n - 1 ) 

and h{a n - 1 ) = 1, i?f(0) = E[X{a n )X{a n )] = a 2TH ' +n , n = 0, 1, . . . , T - 1. Thus 
„ f a(W"'+™ u + n-l<T-2 

^ n (fci + V ; - | a(fe+3 )Tif' + n ti + n— l^T— 1 

Also by straight calculation from (3.18) we have the same result. □ 

Corresponding to the DT-SIM process, {X(a k ),k G Z} with scale I = a T , a > 1, 
T G N there exists a T-dimensional discrete time self-similar Markov process W(t) = 
(W°(t), W 1 ^),. . . , W T ~ l {t)) with parameter space f = {l n ;n G W, I = a T }, where 

W k (l n ) = W k {a nT ) = X{a nT+k ), k = 0,...,T-l. (3.19) 

The elements of the covariance matrix which is defined by (3.6) at points l n and V by (3.10) 
and (3.11) can be written as 

Qf k (l n ,l T ) = E[W j (l n+T )W k (l n )] = a 2nHT E[X{a TT+j )X(a k )} 

= a 2nHT R« (tT + j -k) = a^ HT [h{* T - l )YCf k RnV) (3.20) 

in which Cf k = h(a^- 1 )[h(a k ~ 1 )]~ 1 and R% (.) is defined in (3.10). 

Theorem 3.3 Let {X(a n ),n G W} be a DT-SIM process with the covariance function 
Rn(r). Also let {W(l n ),n G W}, defined in (3.19), be its associated T-dimensional dis- 
crete time self-similar Markov process with covariance function Q H (l n ,l T ). Then 

Q H (l n ,r) = a 2nHT C H RH[h(.a T - 1 )] T , rGW (3.21) 

where h(-) is defined by (3.12) and the matrices Ch and Rh are given by Ch = [CjQj,jfc=o,i,...,T-i; 
where Cf k = h(cp- 1 )[h(a k - 1 )]- 1 , and 



Rh = 



Rg(0) ••• 
flf(0) ••• 



.□ 
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Remark 3.1 It follows from Theorem 3.3 that for each k = 0, 1, . . . ,T — 1 the process 
W k (l n ) = X(a nT+k ) is a self-similar Markov process for n G W. The covariance function of 
the process is 

rf (l n , V) = E[W k (l n+T )W k (l n )] = a 2nHT [h(a T - 1 )] T R%(0), r G W 

where = h(a k - 1 )[h(a k - 1 )}- 1 = 1. 

The introduced T-dimensional self-similar Markov process W(t) with parameter space 
T = {l n ,n G W}, Z = a n , T G N, is the counterpart of the T-dimensional stationary 
Markov process Y(t) = (Y°(t),Y 1 (t),...,Y T - 1 (t)). The spectral density matrix of such 
T-dimensional self-similar process is characterized by the following lemma. 

Lemma 3.4 The spectral density matrix d H (cj) = [dj^(u;)]j ir= o,...,T-i of the T-dimensional 
self-similar process {W(l n ),n G W} ; defined by (3.19), where I = a T has the Markov property 
and is specified by 



2tt 



/i(a'- 1 )(l - e-™ T a- HT h{a T - 1 )) h{a^ l )(l - e-^a HT h- 1 (a T - 1 )) 



where R^(0) is the variance of the process X(-) at point a k and h(a k ) is defined by (3.12). 

Proof: As we mentioned prior to the proof of the Theorem 3.1, we consider Qf r (-,-) at 
discrete points l m and I s where m,seZ, then 

Qf r (l m ,l s ) = E[W j {r +s )W r {l m )] = l 2mH E[W j (l s )W r (l)} = l 2mH Qf r (l s ). 

If the T-dimensional discrete time self-similar Markov process W(-) is sampled at points 
P = a nT , then in (3.7) we have tT instead of r, thus in (3.5) we have nT instead of n and the 
corresponding spectral density matrix of the covariance matrix Q H (l s ) = [Qf r (l s )]j, r =o,i,...,T-i 
is 

d H (u) = [d^( W )] J - r=0) l 1 ... 1 T-l 

where 

d » = t jtl- Hs e— T Qfr(l S ) + ^ £ r Hs e— T Qf r (l s ) ■= df rl H + df r2 H. 

s=0 s=— oo 

Using (3.20), we evaluate d^ rl (co) as 

1 oo 1 oo 

s=0 s=0 



oo 



Mg^WW y (e-^a-^M^- 1 )) 3 . (3.22) 
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Now we verify the convergence of the above summation. By (3.12) we have 



le-^hioF- 1 : 



j=0 j=0 



T-2 

n 

3=0 



E[X(ai +1 )X(ai)} 



Rf(o) 

E[X(a T )X(a T - 1 )} 



By scale invariance of X(-) we have that E[X 2 (a T )] = a 2TH E[X 2 (1)}. Now for j = 0, . . . , 
T - 1 if at least one of the Coii[X(a^ +1 )X(a^)} < 1 then h(a T ~ 1 ) < a TH , and 



^ T a~ HT h{a T - l )\ < 1. 
Therefore the summation on the right side of (3.22) is convergent. Thus the spectral density 



is 



hiai- 1 )^ (0) 



Now we are to evaluate dj^ 2 (aj) as 



1 _ e -iuT a -HT h ( a T-l^ 



1 oo 



s=l 



As Qf r (l~ s ) = E[Wi(l~ s )W r {\)] = l- 2sH E[Wi{l)W r {l s )] = l~ 2sH Q^{l s ), so by a similar 



method, one can easily verify that 



fe(a r - 1 )H:f(0) e iwT a- ifT / l (a T - 1 ) 



2vr/i(^'- 1 ) 1 - e^ T a- HT h{a T ~ 1 ) 
So we arrive at the assertion of the lemma. □ 



Remark 3.2 Lemma 3.4 provides the spectral density of discrete time self-similar Markov 
process {W k (l n ),n G Z} ; defined by (3.19), /or fc = 0, 1, . . . , T - 1 as 



i?f(0)(l-a- 2OT /i 



2HTu2f„,T-l 



a 



L )) 



2tt(1 - 2cos(wT)a- ffT /i(a T - 1 ) + a- 2 ^/i 2 (a r - 1 )) 



Remark 3.3 Using Lemma 3.4, relations (2.16), (2.17) and (3.20), me see that the spectral 
density matrix f(uj) = [fjr(^)]j,r=o,i,...,T-i of a DT-SIM process which is defined by (2.18) is 
fully specified by {Rf {!) , Rf (6) , j = 6,1,..., T-l}. 
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Example 3.2 Here we present the T- dimensional discrete time self-similar Markov process 
corresponding to the simple Brownian motion, described in Example 3.1, as W(l n ) = (W°(l n ), 
W l (l n ), . . . ,W T ~ l (l n )), where W k (l n ) = X(a nT+k ). Now as we mentioned in Lemma 3.5 
we obtain spectral density matrix ofW(l n ). In Example 3.1 we find that h(a^ 1 ) = 1 and 
h{a r - 1 ) = 1 as j,r = 0,1,..., T- 1, h(a T - 1 ) = a TH ' , H' = H - ± and i?f(0) = a 2TH ' +r , 
Rf(ty = oi 2TH +J thus the spectral density matrix ofW(l n ) is 



2TH' 

df» = — 



a a J 



1 - e-^ T a- T / 2 1 - e-^ T a T / 2 
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